Shot noise in ferromagnet-normal metal systems 
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A semiclassical theory of the low frequency shot noise in ferromagnet-normal metal systems is 
formulated. Noncollinear magnetization directions of the ferromagnetic leads, arbitrary junctions 
and the elastic and inelastic scattering regimes are considered. The shot noise is governed by a set of 
mesoscopic parameters that are expressed in terms of the microscopic details of the junctions in the 
circuit. Explicit results in the case of ballistic, tunnel, and diffusive junctions are evaluated. The 
shot noise, the current and the Fano factor are calculated for a double barrier ferromagnet-normal 
metal-ferromagnet system. It is demonstrated that the shot noise can have a nonmonotonic behavior 
as a function of the relative angle between the magnetizations of the ferromagnetic reservoirs. 
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I. INTRODUCTION 

Magnetoelectronic circuits have recently attracted con- 
siderable interest due to their potential for magnetic 
random access memories, sensors, and for fundamental 
studies (Qf-ii^Pjin transport in magnetic and nonmagnetic 
devices. ErBBa A ferromagnetic metal in contact with a 
paramagnet can inject spins inta the paramagnet or de- 
tect spins from the paramagnet .El Spin accumulation and 
its effect on the current-voltage characteristics have been 
studied thoroughly in double barrier ferxomagnet-normal 
metal-ferromagnet (F-N-F) systemsBEl 

An important parameter in view of applications of spin 
injection is the noise to signal ratio. The noise further- 
more provides additional information on the electronic 
structure and the nonequilibrium transport processes and 
that is our main motivation for this study. The calcula- 
tion and experimental detection of shot noise has been 
an actiws subfield of mesoscopic physics during the last 
decade.u Shot noise in systems consisting of normal met- 
als, superconductors, semiconductors and combinations 
thereof combined by tunnel barriers, -diffusive barriers 
and ballistic barriers have been studiedB There have been 
much less attention on the shot noise in hybrid systems 
involving ferromagncts. 

Only a few studies have been carried out on the effects 
of the fluctuations in the spin accumulation on the shot 
noise in ferromagnet - normal metal systems. Bulka et al. 
computed in Ref. ^ the shot noise in F-N-F double barrier 
systems in the Coulomb blockade regime when the mag- 
netizations in the ferromagnetic reservoirs are collincar 
and found pronounced and interesting effects due to the 
interplay of the spin and charge fluctuations . Their re- 
sults in the regime when the source-drain bias is larger 
than the Coulomb charging energy can be understood in 
terms of well-known-results for noipaal-metal systems for 
two spin directions.^ Nowak et alS3 measured electrical 
noise in ferromagnet-insulator-ferromagnet (F-I-F) sys- 
tems in the collinear configurations. They also obtained 
results consistent with a generalization of the results for 
the normal metal systems. Spin-dependent tunneling 
through F-I-F systems in the case of noncollinear mag- 
netization configurations were experimentally studied by 
Moodera et aZ.til with results in good agreement with the 



predictions for the angular dependence of magnetoresis- 
tance by Slonczewski,t3 but no noise measurements were 
performed, ft is our purpose to provide a more general 
discussion of the shot noise in ferromagnet-normal metal 
systems where Coulomb charging effects are negligible. 
The noise power S(ui) is defined as 
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where AJ(t) = lit) — (I(t)) denotes the fluctuation of 
the current I(t) from its average (I(t)). The noise power 
has contributions from I// noise, thermal noise, and shot 
noise. In equilibrium, the thermal noise is given by the 
Johnson-Nyquist formula in terms of the conductance of 
the system. We will consider the nonequilibrium noise 
contributions in the limit when the bias voltage is larger 
than the thermal energy eV ^> k-gT, so that the thermal 
Johnson-Nyquist noise can be disregarded. The 1/ f noise 
dominates at very low frequencies or high bias. 

Although our theoretical framework allows the calcu- 
lation of the shot noise in multiterminal systems, we will 
focus our attention on double barrier systems. Double 
barrier systems have a characteristic time scale corre- 
sponding to the dwell time Td of the electron on the node. 
We will consider the low-frequency shot-noise contribu- 
tion to the noise power when u -C t^ 1 but the frequency 
to is higher than the frequency at which the 1// noise 
dominates. Thus in this regime the low- frequency shot 
noise is the dominant term and we will set u> = in the 
following and disregard the I// noise. 

The Landauer-Biittiker (LB) theory of transport in 
phase-coherent mesoscopic concLuctors gives simple ex- 
pressions for the average currentfij 



and the low-frequency shot-noise power 
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where V is the applied voltage and T n (transmission 
probabilities) are the eigenvalues of tt' composed from 
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the transmission matrix t of the junction at Fermi en- 
ergy. It is assumed that the scattering matrix is en- 
ergy independent in the energy interval determined by 
the chemical potentials of the voltage probes. This is 
easily realized in metallic systems and/or at low bias. 
Of special interest is the Fano factor which is defined as 
the ratio between the shot-noise power and the average 
current: F — S/(2eI). When the transmission of the 
particles of charge e is Poissonian (i.e., random) the shot 
noise is S = 2el, and the Fano factor is F = 1. 

We consider a double barrier F-N-F system as shown 
in Fig. [j]. The system is composed of a chaotic cavity 
connected to two voltage probes at different chemical po- 
tentials through junctions. The left and the right probes 
are ferromagnetic with magnetization directions mx and 
m 2 , respectively. Our theory is applicable if the follow- 
ing assumptions hold: The spin-flip relaxation time on 
the normal metal island is much longer than the typical 
dwell time on the island. The size of each of the two 
junctions is much smaller than the corresponding phase- 
coherence length, but the motion of the electrons inside 
the normal metal island can be elastic or inelastic. The 
electrons trapped in the cavity dwell long enough so that 
the scattering in the cavity leads to a homogeneous and 
isotropic phase-space distribution of electrons inside the 
cavity. This can be achieved, for example, by scattering 
off the bulk disorder (diffusive cavity) or by scattering 
off the irregularities at the surface (ballistic cavity). The 
theory is applicable if the cavity has negligible resistance 
compared to the resistances of the junctions, e.g. the 
voltage only drops across the junctions. Two scattering 
regimes are considered: quasielastic when t c i <§; 7m, and 
inelastic when T- m -C Td, where T m is the inelastic scat- 
tering time, and Td is the dwell time in the cavity. 

We will study the influence of the relative difference 
in the magnetization directions on the shot noise. In 
the case of collinear magnetizations, the results for the 
all-normal metal systems can be easily generalized by 
taking the two spin directions separately into account. 
The regime of noncollinear magnetization configuration 
is more interesting and we will present novel nontrivial 
results for this case. There are, to the best of our knowl- 
edge, no systematic predictions for the Fano factor in this 
regime. We will also present results for the shot noise in 
other transport regimes than previously studied for F-N 
systems. In addition tO|-,th|e tunneling regime, already 
discussed in some limits,EnI9 we will extend the analysis 
to ballistic and diffusive junctions. Different regimes of 
scattering on the normal metal node will also be studied: 
The elastic scattering regime and the inelastic scattering 
regime. 

Let us first review the results for transport through 
a single barrier in all-normal metal systems.El The Fano 
factor in the cases of ballistic (B), diffusive (D), and 
tunnel (T) junctions are 
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In the ballistic and tunnel regimes the result follows di- 
rectly from the LB theory (|^,|3|) , since T n <C 1 for tunnel 
junctions and T n = or T n — 1 for ballistic junctions. 
The 1/3 noise suppression in diffusive junctions is due to 
the bimodal form o£ the distribution function of trans- 
mission eigenvalues. Lj 

Next, let us review the results for all-normal metal dou- 
ble barrier systems in the elastic transport regime. The 



Fano factor is F — (RfFi 
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R2) 3 , where J4 (Ri) is the Fano factor (resistance) of the 
ith junction.li-3 For ballistic, diffusive and tunnel junc- 
tions the Fano factors for double barrier systems are 
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In the ballistic case (f7j), each junction exhibits zero shot 
noise, but the shot noise of the double barrier system is 
finite since the chaotic nodes have a partial (fluctuating) 
occupation of the energy levels in the nonequilibrium sit- 
uation even at zero temperature. Consequently, although 
the direct contribution to the shot noise from a single 
junction vanishes, the effective 'thermal' contribution to 
the noise is finite. The elastic chaotic cavity was recently 
realized as a quantum dot connected to voltage probes 
with two point contactsEj and the Fano factor (Q) was 
experimentally confirmed. For diffusive junctions, the 
Fano factor (|sj) is the same as for a single junction (||) 
because an elastic double diffusive barrier system can be 
seen as a single diffusive junction. For tunnel junctions, 
the Fano factor (||) attains its minimum, F T X 2 = 1/2, 
when the barriers are symmetric. For completely asym- 
metric barriers, the system is governed by the barrier 
with the lowest transmission only and therefore the Fano 
factor is F T l 2 = 1 as in the case of a single tunnel barrier. 

In the inelastic transpoxt regime the Fano factor for 
a double barrier system isEJ F = (i??Fi + R%F 2 )/(Ri + 
R2) 2 ■ Because the energy level occupation numbers in 
the node obey the Fermi-Dirac (FD) distribution with 
a fluctuating chemical potential, the Fano factor of the 
entire system is a linear combination of the Fano fac- 
tors of the two junctions. Ballistic, diffusive, and tunnel 
junctions then yield, respectively, 
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In the case of tunnel junctions, we get the same result 
as in the elastic regime. The Fano factor in the diffusive 
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double barrier system differs from 1/3, which is a univer- 
sal noise suppression factor in elastic diffusive wires. 1/3 
suppression is restored only for a completely asymmet- 
ric system governed by the lowest transmission barrier. 
Finally, ballistic system has vanishing noise because the 
node has a FD distribution of electrons. 

We will generalize the results above to ferromagnet- 
normal metal systems with noncollinear magnetization 
directions. The manuscript is organized in the following 
way. In Sec. || we describe the circuit theory used to 
analyze multiterminal F-N systems and relate the cur- 
rent and its fluctuations through a junction to the ad- 
jacent distribution functions and the scattering matri- 
ces for spin- up and spin-down electrons. We compute 
the Fano factor for double barrier junction s yst ems in 
Sec. [II. Our conclusions can be found in Sec. |v[ 



II. CIRCUIT THEORY 

Nazarov realized that transport in hybrid 
superconductor-normal metal systems can be un- 
derstood in teems of a generalized circuit theory similar 
to Ohm's law£3 In a similar way a circuit theory for 
ferromagnet-normal metal systems was formulated by 
Brataas, Nazarov, and Bauer.Q We extend the circuit 
theory of transport through ferromagnetic-normal metal 
systems to the calculation of shot noise. A (mag- 
neto)electronic circuit is divided into junctions (resistive 
elements), nodes (low impedance interconnectors), and 
reservoirs (voltage probes at thermal equilibrium). 
The theory is applicable when the junctions limit the 
electric current and the nodes are characterized by 
a distribution function which is constant in position 
space (this assumption can be relaxed, see Ref. 17|) and 
isotropic in momentum space. It is assumed that there 
are no spin-flip processes in the nodes. In the following 
discussion the nodes are taken to be normal metal 
chaotic cavities. Each voltage probe is ferromagnetic, 
connected to the nodes through leads carrying current 
and spin current. The current in the circuit dependsppn 
the relative magnetization orientations in the probes B 

The current in the circuit can be_derived by using the 
Keldysh Green function technique!] We obtain the same 
result using LB approach generalized to describe circuits 
with one or more terminals (the nodes) with nonequilib- 
rium spin distribution functions. The assumption of the 
LB theory is that particles exiting a probe have the same 
occupation of the energy levels as the particles inside the 
probe and the occupation numbers of the incoming par- 
ticles are uniquely determined by the probes that sup- 
ply them and the scattering properties of the junctions. 
We generalize the LB expression for the current operator 
from a normal metal node into a junction to a 2 x 2 cur- 
rent operator in spin space in order to also describe spin 
currents collinear and noncollinear to the magnetization 
directions in the adjacent ferromagnets. The 2x2 cur- 
rent operator for the particle flow from node 2 to node 1 
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where a ani (E) creates a spin-a electron with energy E 
leaving the ith node through the nth transverse mode and 
ban i (E) creates a spin-a electron with energy E entering 
the zth node through the nth transverse mode. The cur- 
rent operator is I = ^2 a I aa and the spin-current opera- 
tor is i s = -h/(2e)J2 a p&apI (3a , where cr = {a x ,a y ,a z ) 
is a vector of the Pauli matrices. If the scattering ma- 
trices of the junctions are known then we can express 
annihilation operators b an ,i(E) in terms of the opera- 
tors a an ^(E). Expectation values involving operators 

Q> an i(E) and a an ,i(E) can then be evaluated in terms 
of the distribution functions in the adjacent nodes and 
we can calculate both the current and the shot noiseJlj 
Let us consider a node that is connected to several 
other nodes via junctions. In the low frequency, long- 
time limit u> <C 1 /Td the charge and spin in the node are 
conserved at every instance of time 



(14) 



where 7j is the current flowing out of the node through 
the ith junction at time t and the index i runs over all 
junctions connected to the node. Due to the charge and 
spin conservation the distribution function on the 

normal metal node has two contributions: A stationary 
part f^(E) and a small fluctuating part 5/m(E, t) caused 
by the discreteness of the charges and spins that enter 
and leave the island. Using the isotropy assumption in 
momentum space, the stationary part f^(E) is defined 
by (at m (E)a Pn (E')) = 6 nm S(E - E')f^ a (E), where <J am 
creates a spin-a particle in the mth quantum state of 
the node. The distribution function /n{E) is a 2 x 2 
Hermitian matrix in spin space to allow a nonequilibrium 
spin accumulation on the normal meJal_node. 

Following Beenakker and Buttikerj3'til we separate the 
time-dependent current fluctuations A/j(t) in the ith 
junction into two contributions. The first contribution 
8Ii(t) is due to the intrinsic noise in the junction if the 
fluctuating term 8fpf(E, t) of the distribution function on 
the node is neglected. This is the only relevant term for 
the shot noise for a single junction connected to reser- 
voirs that have stationary distribution functions. The 
second contribution is due to the small fluctuation in the 
occupation of the energy levels on the node 8f]sr(E, t). 

Spin and charge conservation ( p"4] ) dictate the average 
occupation of the energy levels in the node in terms of 
the average distribution functions in the adjacent nodes. 
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From Eq. ( |l4| ) the total fluctuations of the current vanish 
]TAi 4 (*) = 0. (15) 

i 

Consequently, small fluctuations of the occupation of the 
energy levels are needed to compcnsate-for the intrinsic 
current fluctuations in the junctionsJ§E-J The total fluc- 
tuation AIi(t) is a sum of these two contributions: 
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where j runs over the nodes adjacent to the ith junc- 
tion. The calculation of the term proportional to Sfj(t) 
resembles the case of all-normal metal systems, but is 
now generalized to describe the 2x2 spin-space matri- 
ces. The calculation of this term is straightforward, but 
tedious since it requires solving for the 2x2 distribution 
functions in the circuit. The explicit result for the total 
fluctuation in the current A Ii (t) in terms of the intrinsic 
fluctuation in the current through the separate junctions 
8Ii(t) is shown below in Eq. ( |30| ) f° r symmetric double 
barrier systems. 

The ferromagnetic reservoirs are in thermal equilib- 
rium characterized by equilibrium FD distribution func- 
tions independent of time. In contrast, the normal metal 
islands are isolated from the thermal baths and have 
(fluctuating) nonequilibrium distribution functions re- 
quired for current conservation. The normal metal nodes 
can be viewed as fictitious probes £-3 They are character- 
ized by an isotropic and homogeneous distribution func- 
tion like the true probes (i.e., probes in local thermal 
equilibrium), but the (fluctuating) distribution function 
does not have to be of the FD form. 



A. Average current 

Naturally, we obtain here the same result for the av- 
erage current by using the LB formalism as from the 
Keldysh formalisms since they are identical when there is 
no inelastic scattering in the junctions. For completeness 
we show the result for the average current (per unit of 
energy, at a given energy) between a ferromagnetic jjiode 
and a normal metal node at the normal metal sidefl 

3 T u t (/f - /W)w T + g l u l (fF - /n)u 1 



-g^vJ Jnu^ — fNvJ 
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where u' = (1 + a ■ m) /2 and w- = (1 — a • m) /2 are pro- 
jection matrices, /jv is the average distribution function 
in the normal metal node, fp is the equilibrium distri- 
bution function in the ferromagnetic probe, and the con- 
ductances are defined in terms of the reflection matrices 
for electrons incoming from the normal metal node 
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For energies much less than the Fermi energy, the conduc- 
tance (|l8|) is energy independent. Therefore, the elastic 
and inelastic transport regimes are equivalent for the av- 
erage current. However, the shot noise differs in these 
transport regimes, since it is sensitive to the the energy- 
resolved distribution.^ 



B. Shot noise 

The starting point in calculating the noise is the rela- 
tion between the fluctuation of the current in a junction 
and of the distribution function in the fictitious probes 
adjacent to the junction (Jlq ). The first term on the 
right-hand side of the Eq. (|16|) is a spontaneous Langevin 
source. The second term ( P_6| ) can be found in terms of 
the Langevin source by using Eq. (15h for the instanta- 
neous conservation of spin and Eq. (|l7p for the current. 
The problem of calculating the noise is then reduced to 
finding the correlator of the Langevin sources Sli in terms 
of the scattering matrices of the junctions and the distri- 
bution functions in the adjacent nodes and reservoirs. 

The correlators of the Langevin sources can be found 
by generalizing the LB theory for all-normal metal sys- 
tems. We expand the current in terms of the creation 
and annihilation operators, using Eq. (|l3|), and then de- 
termine particle and spin fluctuations 



- (<4(£i) 0/J i(3»)) (a\ m (E 3 )a Sn (E 4 )) 
= 5 kn S lm S(E 1 - E 4 )S(E 2 - E 3 ) 
xf s n a (E 1 )\l-f^{E 2 )] , (19) 



where the Greek subscripts denote spin indices, and the 
Latin subscripts include both the transverse mode index 
and the probe label. However, the distribution function 
in the probe /„ is independent of the transverse mode 
and n here denotes only the probe label. 

A similar analysis was used to calculate shot noise in 
nonmagnetic circuitsJi-3 Our treatment of spin transport 
is a generalization to include the spin degree of freedom. 
In addition to instantaneous conservation of charge in the 
nodes, we use instantaneous conservation of spin. Both 
conditions are satisfied when the frequency is lower than 
the inverse dwell time of a particle in the node, lot& <C 1. 

When the inelastic scattering time is much longer than 
the transport dwell time, the energy of a quasiparticle is 
conserved. In this regime both current and noise can be 
calculated at each energy level and then integrated over 
the energy to get the final result. 

The inelastic regime is achieved when the inelastic 
scattering time is shorter than the transport dwell time. 
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When a particle is trapped in a node, it is assumed to 
drop to the lowest energy state allowed without a spin 
flip. As a result, there will be a direction of spin accu- 
mulation in the node with FD distribution for electrons 
polarized parallel and antiparallcl to it (with different 
chemical potentials fv and fj,^, respectively). The direc- 
tion of spin accumulation and its magnitude are deter- 
mined by the transport rates integrated over all energies. 

C. Junction conductance and shot-noise matrices 

Because the current ( |l3[ ) is expressed as a linear com- 
bination of products of two creation and annihilation op- 
erators, the current in the circuit is completely deter- 



mined by 2 x 2 conductance matrices (|18|), elements of 
which can be expressed as traces of two reflection matri- 
ces. Shot noise is quadratic in current, and we expect 
that it will be governed by similar traces with four re- 
flection and transmission matrices. We show that this is 
indeed true with some simplifications. First, because fer- 
romagnetic probes are in local thermal equilibrium, the 
new independent parameters can be expressed in terms of 
the traces of four reflection matrices only. Second, only 
the traces of the form Tr [r a (r a yr"(r" )'] enter the shot 
noise, where a, a' , 0, and 0' denote the spin. The new 
set of parameters describing the shot noise of a junction 
can therefore be grouped into a 4 x 4 Hermitian matrix 
s: 
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(20) 



where the trace is taken for each element of the 4x4 
matrix (each of which is an M x M matrix in the basis of 
the M transverse modes) inside the square brackets, i.e., 
the result of the operation is a 4 x 4 matrix of complex 
numbers. 

Only four off-diagonal elements and three diagonal ele- 
ments are independent in the shot-noise matrix s. There- 
fore, in general, we introduce 7 new parameters: The real 
parameters Sj, s-fj,, sj,, and the complex quantities s+, 
S-, so, s , 



( «T s + s + s o 
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These parameters ( f2l| ) and the conductance parameters 
( |l8|) completely determine the current and noise in the 
system. Scaling of all r and s matrices of a given system 
by the same factor does not change the Fano factor. 

We proceed with the explicit evaluation of the con- 
ductance g and shot-noise s matrices in the cases of tun- 
nel, diffusive and ballistic junctions below and then apply 
our theory to two-terminal double barrier ferromagnetic- 
normal-ferromagnetic (F-N-F) systems. The results for 
the conductance matrix g have been derived in Ref . [3] and 
we briefly reiterate these results here for completeness. 
The shot-noise matrix s has not been studied before. 



Ballistic junctions 



Because the transverse momentum is conserved in bal- 
listic junctions, the reflection matrices (r^ and H - ) are di- 
agonal in theJiasis of the transverse modes. Simplifying 
the situation,E£l we assume that the diagonal components 
of the reflection matrices can attain only two values: Full 
transmission or no transmission 1. From Eqs. (||), (||) 
we see that the Fano factor for a single ballistic junction 
vanishes. It is not the case in a double barrier system, as 
discussed in the Introduction. 

The conductance g and shot-noise s matrices for a sin- 
gle junction are in this case 




g/g 
s/g 



where g = [g^ + g^)/2 is the average conductance and 
p = (g' — g^)/(g* + g^) is the relative polarization. If 
all the junctions are the same, noise suppression will de- 
pend only on the relative polarization p and the relative 
magnetization orientations. 

2. Diffusive junctions 

The conductance for diffusive junctions can be found 
with different methods. One possibility is to solve the 
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diffusion equation with the proper boundary conditions 
as was done in Ref. || to find the mixing conductance 
g*+. An alternative approach is to evaluate the conduc- 
tance matrix by using random matri x theory (RMT) in 
the semiclassical approximation .I13E3 We will use the lat- 
ter approach here, because it considerably simplifies the 
calculation of the shot-noise parameters. 

We assume that the junction consists of two connected 
parts as shown in Fig. ^|. The first, normal metal part 
is described by a single scattering matrix for both spin- 
t and spin- J. electrons. The second, ferromagnetic part 
is described by two independent scattering matrices, one 
for spin-f and one for spin- j electrons. It is assumed that 
there are no correlations between the scattering matrices 
of the spin-1" and spin-j electrons in the ferromagnetic 
part. However the up- and down-spin parts of the to- 
tal scattering matrix of the combined normal metal and 
ferromagnetic system are correlated since both spin di- 
rections see the same scattering centers in the normal 
metal part. Scattering at the F-N boundary is disre- 
garded since it is assumed that the total resistance is 
dominated by the diffuse normal metal and ferromag- 
netic metal parts of the junction. The total reflection 
matrix r a for spin-a electrons can then be found by con- 
catenating the normal metal and ferromagnetic parts 



rN + t' N r%Y(r' N rp n 



n=0 



Tt N 



(22) 



We apply the standard polar decomposition for re- 
flection matriceSr-ef the spin-a electrons inside the fer- 
romagnetic part:B r F = iU a ^/TfU' a , where U a , U' a 
are M x M unitary matrices, which in the isotropic 
approximation 2 .!! are uniformly distributed in the group 
U(M). In the case of time-reversal symmetry, these ma- 
trices are related by transposition, U' a — U aT . T F is 
a diagonal matrix containing the eigenvalues of t F (t F )^ . 
After inserting these reflection matrices into the expan- 
sion (p2|), we perform the averaging of traces entering 
conductance g and shot-noise s matrices over the scat- 
tering matrices in the ferromagnetic part. To this end 
we use-the semiclassical result for traces of matrices from 

rmtJi3 

In order to perform the averaging, expressions of the 
following form must be evaluated. 

p = Tr[(AaBf3 ■ ■ ■ r]C)(F^J ■ ■ • attftypt)] ; (23) 

where A, B, C, D, E, F are some fixed M x M matrices 
and a, (3, 77, 7, 5, u> are unitary-jnatrices uniformly dis- 
tributed in U(M). It was showr£3 that if the ordered sets 
a, P, and 7, S, . . ., u> are such that any two closest 

neighbors are independent (or related by transposition, 
i.e., a = f3 T ), then, to the leading order in M, this trace 
is nonzero only if A = D, B = E, . . ., C = F and a = 7, 
(3 = S, . . ., r) = lo. If these conditions are satisfied, the 
trace ( p3| ) equals (to the leading semiclassical order in 
the inverse number of transverse modes 1/M) 

p = M±-Tr[AAl]±-Tr[BBi] ■ ■ ■ ^[C&] . (24) 



This is sufficient to calculate both the conductance g 
and the shot-noise s matrices in terms of the conductance 
of the normal metal part and up- and down-spin conduc- 
tances in the ferromagnetic part, after we resolve two 
problems. First, we encounter the task of calculating the 
trace Txffjvijyijvtjy]. Second, the two s matrix elements 
Tr [l_ r T( r T)t r T( r T)t] and^ = Tr[l- 



pose a challenge, because some of the terms in their 
expansions [following Eq. (53)] contain correlations be- 
tween random matrices which invalidate the assumptions 
of the semiclassical result ( [m]) . We calculate these three 
quantities using the Boltzmann-Langevin approach. It 
was shown in Ref. [2^ that in all-normal metal systems 
the Fano factor has a universal 1/3 suppression in diffu- 
sive junctions of arbitrary shape, dimension, and con- 
ductivity distribution. This is directly applicable to 
the case at hand because in evaluating Tr^Tv^jv^^JV 1 



and Tr[r Q (r Q )tr Q (r 
Therefore we find 



a y] we deal with collinear transport. 



Tr[tTft' N t N t' N ] 

h 



3 9N - 



(25) 
(26) 
(27) 



Now we use the expansion (E2|), the semiclassical re- 
sult (||), and the results @, (§§), and @ to find the 
conductance and shot-noise matrices 



9/9 



3 ~/ 



9n/9 
/2(l+p) 



l+P 9N/g 
l-p 

l+P 1 
1 



(28) 



V 



V 
P 



1 







-9N/g 



-9/9N 





1 - 


fp 1 + 


P 2(1 




1 


1 2X 









1 




(1 





1 






1 


1 0/ 

















1 


-p 2 1 


-p 2 





1 


-p 2 1 


-p 2 


Vo 











p 
p 

-p) 



(29) 



M 



where gjy is the conductance of the normal metal part of 
the junction, g = (,g T +.g i )/2, andp = (g 1 -g l )/{g^ +g l ). 
The spin-dependent conductance g^ = {g\r9N) / {9p+ 9n) 
and g F = (g F gN)/(g F + ffjv) are given by Ohm's law in 
terms of the spin-dependent conductance of the ferro- 
magnet (g F and g F ) and the conductance of the normal 
metal (gjv)- Two parameters govern the Fano factor if 
all junctions are the same, the relative polarization p and 
the fraction of the conductance of the normal metal part 
to the average conductance gj\r/g- The result for the 
conductance matrix (|2^) using RMT agrees with the cal- 
culation using the diffusion equation in Ref. ^. 



7 



3. Tunnel junctions 

The transmission coefficients are exponentially small 
in tunnel junctions. We expand f and s in terms of the 



small quantities 5r° 



and keep only the 



lowest order nonvanishing terms. We also assume that 
the reflection coefficients have random phases. First- 
principles bancUsitructure calculations confirm this for re- 
alistic systems.Ej An important result can be drawn from 
the randomness of the phases of the reflection coefficients: 
The imaginary part of the mixing conductances and, 
similarly, the shot-noise parameters vanishes. Using this, 
we can express the conductance and shot-noise parame- 
ters as 



», ( 1+P 1 
9/9 = 11- 



P 



s/g = 



l + 2pl+pl+p 1 
1 +p 1 1 1 -p 
1+p 1 1 1 -p 
1 l-pl-pl-2p 



The parameters only depend on the average conductance 
g and the relative polarization of the junctions p. 



III. SYMMETRIC F-N-F DOUBLE BARRIER 

The theory developed above can be used to ana- 
lyze the shot noise in complicated many-terminal de- 
vices containing many normal metal nodes and ferro- 
magnetic reservoirs. Such an analysis has already been 
performed for the current through two-terminal double 
barrier ferromagnet-normal metal-ferromagnet systems, 
a novel three-terminal ferromagnet-normal metal- fer- 
romagnet "spin-flip" transistor and a generalization of 
Johnson's three-terminal spin-transistor.tl We have seen 
that, in general, there are 4 conductance parameters and 
11 shot-noise parameters that completely describe each 
junction in the system. The Fano factor defined as the 
ratio between the shot noise and the current through a 
junction thus depends on 15N — 1 parameters, where N is 
the number of different junctions. (The Fano factor is in- 
variant under the scaling of all parameters and, therefore, 
we can set any one of the 15iV parameters to be unity 
and correspondingly scale the remaining 15-/V— 1 param- 
eters.) Clearly, we cannot explore all possible systems, 
but we will here illustrate the usefulness of the semiclas- 
sical theory of shot noise to study the most simple sys- 
tem of spin accumulation: A two-terminal double barrier 
ferromagnet-normal metal-ferromagnet system. We will 
furthermore assume that the system is symmetric and 
study the regime of ballistic, diffusive, and tunnel junc- 
tions. Besides, we will investigate the elastic and inelastic 
transport regimes separately. The number of indepen- 
dent parameters (defining the angular dependence of the 
Fano factor) then reduces to 1 (ballistic), 2 (diffusive), 
and 1 (tunnel), which considerably simplifies the analy- 
sis. The results for the current and the shot noise in the 



collinear configuration, when the magnetizations of the 
left and the right probes are parallel or antiparallel, can 
be easily deduced from a two channel model of the corre- 
sponding all-normal metal system. All our results in the 
elastic and inelastic transport regimes and for ballistic, 
diffusive, and tunnel junctions agree with these known 
results. The results when the magnetizations are non- 
collinear are novel. We demonstrate that the Fano factor 
can exhibit a nonmonotonic dependence on the relative 
angle between the magnetizations in the ferromagnets. 
This can be used experimentally to obtain more infor- 
mation about the spin accumulation and the nature of 
the junctions including their polarizations. 

We assume that both ferromagnetic reservoirs are held 
at zero temperature but at different chemical potentials 
Ho and + \eV\, where V is the applied constant volt- 
age bias. We denote the energy-averaged distribution 
function in the normal metal island for energies between 
the chemical potentials of the reservoirs, /jy is a function 
of the conductance matrix g, the relative angle between 
the magnetizations of the ferromagnets, and the applied 
voltage bias and can be found by the average charge and 
spin conservation in the island as outlined in Ref. || and 
is the same in both clastic and inelastic regimes. 

We use Eq. ( [i"7j ) for the current flowing through the left 
and right junctions and Eqs. (|l6|), ( filf ) for the spin and 
charge conservation to obtain the form of the fluctuation 
of the node distribution function <5/at. This fluctuation 
leads to the total fluctuation in the current given by 



with 

Vn -- 



AI 1 {t) = Sp[V 1 6h{t) + V25i 2 {t)} 



; , V (m 2 - mi)r/ r + (m 2 x nnjfjj 



(30) 



■1 



|fy| 2 cos 2 - 



2 , r\ r sin 2 

where r\ r (rji) is the real (imaginary) part of the relative 
mixing conductance rj — 2g^/(g^ + g^). V n is related 
to the average zero temperature distribution function /jy 

by V n = (A - h)(fn - fx), where /„ is the distribution 
function in the nth ferromagnetic reservoir for energies 
between the chemical potentials of the reservoirs. (Be- 
cause the reservoirs are held at thermal equilibrium at 
zero temperature, the allowed values for /„ are 1 and 
0-the unit 2x2 matrix and the zero 2x2 matrix) Sp 
denotes the trace in spin indices to distinguish it from 
the trace in the space of transverse channels which we 
denoted Tr in the preceding discussion. Eq. (^) is valid 
in both the elastic and the inelastic regimes when the 
conductance matrix parameters are energy independent 
on the scale defined by the voltage bias. 

Using Eq. ( p0| ) we express the time-dependent cur- 
rent fluctuations in terms of the Langevin sources 5h (t) . 
Then we apply Eq. (19) to find correlators for these 
Langevin sources and finally find the shot noise from Eq. 

(0): 



o„2 r)io + \eV\ 

S = — dESp 



S 2 (E) 
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S n (E) = Y<9 a0 {fN(E){l - f N {E)){2V n uiV n ul 

-uiV n V n v%) - [f N (E) ~ f n fvP n V n V n u a n } 

aa'Pp' 

XuiV n ui[f N {E) - fn]vgV n } . 



(31) 



In each sum a, a', /3, and (3' are spin indices, u" = 
(1 + am n ■ <x)/2 are projection matrices corresponding 
to the magnetization direction m„ of the nth reservoir, 
fw(E) is the distribution function on the normal metal 
island and g a ^ = Tr[i - r a (r^] and s°' 5a 'P' = Tr[i - 



(r@ Y] are elements of the conductance matrix 



g (|18|) and shot noise s matrix fl20|). 

All the terms in the sums of Eq. (^Tj) are not indepen- 
dent: The elements of the first sum can be grouped to 
form a 2 x 2 Hermitian matrix while the elements of the 
second sum can be grouped to form a 4 x 4 Hermitian 
matrix with the same structure as the shot-noise s matrix 
©. 

In the following subsections we directly apply Eq. ( |31| ) 
to the cases of ballistic, diffusive, and tunnel junctions in 
both elastic and inelastic regimes. 



A. Elastic transport 

1. Ballistic junctions 

The Fano factor for the collincar configurations of the 
ferromagnetic magnetizations can be found by applying 
the result for all-normal metal systems (^) for two inde- 
pendent spin channels 



F{9 = 0) 
F(9 = tt) 



1 

4' 

l-p 2 



The Fano factor is decreased in the antiparallel configu- 
ration 9 — tt when each spin component experiences an 
asymmetric double barrier system. The angular depen- 
dence of the Fano factor is shown in Fig. || for a set of 
different polarizations p. 

For polarizations p below the critical value p c = 1/3, 
the angular dependence of the Fano factor is monotonic. 
When p exceeds p c , there is a maximum in F which con- 
tinuously (and monotonically) increases from 9 = to 
9 = tt as p increases from p c to 1. F max (p) monotoni- 
cally increases from 1 /4 to 1 when p is increased from p c 
to 1. The position of the maximum thus gives another 
independent measurement of the polarization p together 
with the ratio F(0 = %)/F(6 = 0). 

It is interesting to note the limiting behavior 
of F when 9 approaches tt and p approaches 1: 
linip^i lime^Tr F{9,p) — and lim^^ lim p _>i F(9,p) = 



1. Thus, for half- metallic ferromagnets (p — 1), there is 
a sharp drop in F from 1 to when 9 approaches tt. 

In the exact antiparallel configuration, 9 — tt, the 
Fano factor vanishes when the system is close to being 
half- metallic (p is close 1) since the transport properties 
are governed by the barrier with the lower transmission 
which has a vanishingly small Fano factor for ballistic 
junctions, p « 1 means that one barrier has vanish- 
ing conductance and, therefore, the current and the shot 
noise vanish in the system. Relaxing the assumption 
of antiparallel magnetization when the relative angle is 
slightly below tt, the Fano factor increases to 1 as some 
of the low transmission channels form. These low trans- 
mission channels determine the transport properties and 
the system effectively behaves as a single tunnel barrier 
system with F w 1. 



2. Diffusive junctions 

As discussed in the Introduction, the Fano factor has 
the universal value of 1/3 in the collinear configuration 
for diffusive junctions. Thus we have from Eq. (^) for 
two spin-channels 

F(9 = 0) = 
F{9 = tt) = I. 

The 1/3 suppression also holds for arbitrary angles when 
p = so that there is no spin accumulation. In gen- 
eral, there is an angular dependence of the Fano factor 
since the spins are coherent on the normal metal por- 
tion of the junctions. We show this dependence in Fig. 
H for gN/g — 4 and several values of p. When p is fi- 
nite, there is a minimum in F(9). The position of the 
minimum is relatively insensitive to the value of the po- 
larization p. When gN/g increases, this minimum shifts 
to larger angles and the Fano factor flattens in the limit 
of gN I g oo (in this limit the junctions are fully ferro- 
magnetic): F(9) — 1/3 for 9 < tt with a singularity at 
9 = tt when F has a sharp minimum which deepens to 
zero for p — 1. This mimics the behavior of the ballistic 
systems. The Fano factor identically equals 1/3 in a spe- 
cial case of p = 1 and gN/g — 2 when gp/gN oo and 
9f/9n — ¥ 0, where gp is spin-a conductance of the ferro- 
magnetic portion of the junction [for a diffusive junction, 
in general, g a = g F gN/(g F + gw)]- 



3. Tunnel junctions 

In this regime, the Fano factor displays a simple ana- 
lytic dependence on the angle 9 for any polarization p: 



m = 5(1 



2 ■ 2 17 

p sin — 
y 2 



(32) 
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The angular dependence is thus a simple monotonic in- 
terpolation between the values of the Fano factor in 
the collinear configurations [from Eq. (^) for two spin- 
channels] 



For completeness, we plot the corresponding Fano factor 
in Fig. ||. 

B. Inelastic transport 

1. Ballistic junctions 

The collinear configuration in this regime gives a van- 
ishing Fano factor 

F(e = o) = o, 

F(9 = tt) = 0. 

Also when the polarization p = 0, the Fano factor F van- 
ishes for any angle 9. But in the case of noncollincar 
configuration and nonzero polarization p, the Fano fac- 
tor F is always finite. This can be understood as follows: 
In general, the time-dependent distribution function of 
the normal metal node has a FD form for the two spin 
orientations along some (time-dependent) direction u(t). 
When there is a spin accumulation on the node, the chem- 
ical potentials corresponding to the electrons with spin 
parallel and antiparallel to u(t) are different. The elec- 
trons are injected into the normal metal island with spin 
orientations along the magnetization direction of the cor- 
responding reservoir, mi or rri2 . This means that these 
projected components of the distributions functions may 
be partially occupied. Partial occupation of the states 
effectively results in a 'thermal noise' contribution to the 
shot noise of the system. Consequently, the Fano factor is 
finite for intermediate angles 9. This mimics the mecha- 
nism for nonzero shot noise in the all-normal metal elastic 
double barrier system with ballistic junctions which we 
discussed in the Introduction. 

The analytic form of the Fano factor F(9) is relatively 
simple: 

= p 2 sing cos § [1+Jpj (1 + 2 sinf)] 
8(l-b|sin 2 f)(l + Hcos 2 |) 2 ' 

In Fig. |^ we show the angular dependence of the Fano 
factor for a set of different polarizations p. When p in- 
creases from to 1, F(9) yields a single maximum which 
monotonically shifts from 6 = to 6 = iz. -F max (p) mono- 
tonically increases from to 1, when p is varied from 
to 1 . The limiting behavior found in the elastic regime is 
also valid in the inelastic regime. 



2. Diffusive junctions 
In the collinear configuration the Fano factor is 

F(9 = 0) = i, 

F(9 = n) = I+J-. 

The 1 /6 suppression also holds for any 9, when p = 0. We 
show the angular dependence of the Fano factor for a set 
of different polarizations p and gN/g = 4 on Fig. ^. For 
nonzero p, F{9) has two local extrema, one minimum and 
one maximum, the minimum corresponding to a higher 9. 
As in the elastic case, the angles corresponding to the ex- 
trema in F{9) are relatively insensitive to the value of the 
polarization p. The function F(9) becomes monotonic for 
sufficiently small gN/g- In the opposite limit, when this 
ratio goes to infinity, gN/g — * oo, the junctions become 
completely ferromagnetic. In this limit, F(9) = 1/6 for 
9 < 7r with a jump to (1 + p 2 )/6 at 9 = it much like in 
the elastic regime: Before the steep rise, the Fano factor 
F has a sharp minimum at 9 = tt which deepens to zero 
for p = 1 . 

3. Tunnel junctions 

As was noted in the Introduction for the case of all- 
normal metal systems, the elastic and inelastic transport 
regimes yield the same Fano factor. We obtain that this 
remains valid for F-N-F systems even in the case of non- 
collincar configuration, i.e. F is given by Eq. (p^). 

IV. CONCLUSIONS AND DISCUSSIONS 

A semiclassical circuit theory for the shot noise in mul- 
titerminal ferromagnet-normal metal systems has been 
developed. The current is governed by 4 real parameters 
for each contact and the shot noise is governed by 11 
real parameters. The shot-noise parameters have been 
evaluated for diffusive, tunnel, and ballistic junctions. 

The circuit theory has been applied to symmetric dou- 
ble barrier F-N-F systems. Both the elastic and inelastic 
regimes have been studied. In each of the regimes, a 
complex nonmonotonic angular dependence of the Fano 
factor was found in the case of diffusive junctions. In 
the case of ballistic junctions with large conductance po- 
larization, a sharp drop of the Fano factor from 1 to 
was found when the magnetization of the ferromagnets 
approaches the antiparallel configuration. Simple sinu- 
soidal angular dependence of the Fano factor was found 
for the tunnel junctions in both the elastic and inelastic 
regimes. 

We have in this work disregarded spin-flip processes. 
Spin-flip scattering can occur both in the junctions and 
on the nodes and these different processes will influence 
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the current and the shot noise in different ways. The 
spin-flip processes on the nodes can be disregarded when 
the spin-flip relaxation time is longer than the trans- 
port dwell time r B f 3> r^. For shorter spin- flip scatter- 
ing times, the angular dependence of the average current 
is governed by a reduced effective polarization and a re- 
duced mixing conductance.!?! We expect similarly that the 
shot noise is governed by a reduced effective polarization 
and a reduced effective mixing conductance. The antici- 
pated effect of spin-flip scattering in the normal metal 
node on the Fano factor is thus mostly quantitative; 
the effective polarization and mixing conductance are re- 
duced making the angular dependence weaker, and the 
angular dependence eventually vanishes when r s f <C T<j. 
(A discussion of the spin-flip scattering time r s f vs the 
transport dwell time Td has been given by Brataas et al. 
in Ref. [?]). Spin- flip scattering in the junctions can be 
disregarded when the junction is smaller in the transport 
direction than the spin-flip scattering length. If that is 
not the case, spin-flip scattering in the junction can have 
a qualitative effect as well as a quantitative effect both 
on the current and the shot noise. First, we also expect 
that the effective polarization and mixing conductance 
decrease with increasing spin-flip scattering in the con- 
tact. There is also a second, qualitative effect: In the 
absence of spin-flip scattering in the junctions, the cur- 



rent between a ferromagnetic node and a normal metal 
is determined by a 2 x 2 conductance matrix ( |l8|) and 
the fluctuations in the current are governed by a 4 x 4 
shot-noise matrix (pfj|). However, these matrices do not 
completely determine the current and shot noise in the 
presence of spin-flip scattering in the junction since more 
combinations of the scattering matrices become accessi- 
ble (for the calculation of the average current see e.g. the 
appendix in Ref. ||) . Consequently, an expanded basis of 
conductance and shot-noise parameters must be used and 
the angular dependence of both the current and the shot 
noise can be qualitatively different than in the case of no 
spin-flip scattering in the junctions. 
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FIG. 1: The two-terminal device consisting of a normal metal node (AT) attached to two ferromagnetic reservoirs (Fl and F2) 
with arbitrary relative magnetization direction 9. A source-drain bias V is applied between the ferromagnetic reservoirs and 
time-dependent currents Ii(t) and l2(t) flow into the normal metal island from the reservoirs Fl and F2, respectively. The 
contact between the ferromagnetic Fl (F2) and the normal metal node is characterized by the conductance matrices gi (32). 
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FIG. 2: A contact between a ferromagnetic node and a normal metal node. The current is evaluated at the normal metal side 
(dotted line). The transmission coefficient from the ferromagnet to the normal metal is t' and the reflection matrix from the 
normal metal to the normal metal is r. 
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FIG. 4: Angular dependence of the Fano factor F in F-N-F systems in the elastic regime with diffusive junctions. The results 
are shown for different polarizations p = (g T — / (<? T + g^. I g = 4. 
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FIG. 5: Angular dependence of the Fano factor F in F-N-F systems with tunnel junctions. The results are shown for different 
polarizations p = (g T — g^ / (g T + . 
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FIG. 6: Angular dependence of the Fano factor F in F-N-F systems in the inelastic regime with ballistic junctions. The results 
are shown for different polarizations p = (g T — / (g T + g^. 




FIG. 7: Angular dependence of the Fano factor F in F-N-F systems in the inelastic regime with diffusive junctions. The results 
are shown for different polarizations p = (g T — / (g T + g^. I g = 4. 



